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ABSTRACT 

Reconnection shocks in a magneticaUy dominated plasma must be compressive. Non- 
thermal ion acceleration can occur across built-in slow shocks, and across outflow fast 
shocks when the outflow is supermagnetosonic and the fleld is line-tied. Electron acceler- 
ation may be initiated by injection from the dissipation region. Reconnection and shock 
acceleration thus cooperate and non-thermal acceleration should be a characteristic fea- 
ture. 
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Magnetic reconnection converts magnetic energy into particle energy as oppositely 
magnetized flows merge across a thin dissipation region (DR). An upper limit to the (2-D) 
reconnection rate comes from the Petschek model (PK) originally proposed to explain rapid 
solar flare energy release [1]. In PK, the DR length << held gradient length and magnetic 
tension in the outflow thrusts the plasma from the DR. Slow shocks are built explicitly 
into PK, across which the plasma flow changes abrubtly from inflow to outflow. The PK 
tension thrust provides faster reconnection than the Sweet-Parker (SP) type models for 
which the DR length ~ fleld gradient length [2] . 

Though slow shocks are built into PK, similar structures are seen numerically even 
when the DR >> than that of PK [3]. In addition, fast shocks can be present when the 
outflow is supermagnetosonic and the outflow boundary fleld is kept flxed [4]. Although 
real astrophysical plasmas are compressive, most work has employed the incompressible 
approximation to facilitate a global solution. Compressible reconnection models have also 
focused on perturbative solutions and the rate [5], not on the particle spectrum. Here we 
consider the effect of compressive shocks near a DR on the spectrum. 

Slow shocks, unlike fast shocks, have not been extensively simulated, and their po- 
tential for power law acceleration has not been addressed as an important outcome of 
reconnection. The shocks may accelerate ions directly from a thermal distribution and 
the DR may provide the injection electrons required for non-thermal electron acceleration. 
Since much of the flow of a reconnecting region passes through the shocks, reconnection 
may possibly explain sustained non-thermal features in magnetically dominated astrophys- 
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ical phenomena. 

We first solve the jump conditions across the built-in slow shock for the compression 
ratio, given a magnetically dominated inflow. We find the parameter space for which the 
outfiow must be supermagnetosonic and relate the the compression ratio across the slow 
shock, Ts, to that across the fast one, r/. We then discuss the implications for shock 
acceleration. 

The non-relativistic MHD jump conditions for mass, momentum, and energy are [6] 



PlVln = P2V2n, (1) 

Pl^ln + ^1 + BI/StT = P2^in + + BI/StT, (2) 

Plf InVlt - -Bi„Bit/47r = P2V2nV2t " B2n^2t/4:TT, (3) 

{l/2)p^vlvin + r(r - ly^PiVin + {Bl/47r)vin - Vi • BiSi^/47r = {l/2)p2vlv2n 

+r(r - l)-^P2V2n + {Bl/ATT)v2n - ^2 ' ^2B2n/^T^, (4) 



where B is the magnetic field, v is the velocity, P is the pressure, p is the density and F 
is the adiabatic index. The subscript 1(2) refers to the up (down) stream region, and the 
subscript n{t) refers to the normal (tangential) components. The electromagnetic jump 
conditions for an ideal plasma are given by 

Bin — B2n (5) 



(vixBi) = (V2XB2). 

3 



(6) 



The shock is ± to the n, y plane as shown in Fig. 1. We assume the switch-off condition, 
B2y — 0, and also that vi/|f i| • y << 1 [8]. Define c = cos9, § = sin9 and i = tan9 where 
$ is the angle between the downstream flow and the shock normal. Define ci = cos(f)i, 
si = sin(f)i and ti = tan(f)i where 0i is the angle between the upstream field and the shock 
normal. The configuration of Fig. 1 is then described by 

Vln = -Vl, Bin = -BiCi, Biy = -BiSi, V2n = -V2C, V2y = V2S, B2n = -B2, Viy = B2y = 0.| 

(7) 

For r = 5/3 and f3\ = afg/v^^ << 1, where au and via are the inflow sound and Alfven 
speed, plugging (7) into (l)-(6) gives 

t? = 2(r,-l)(r,-4)/(5r,-2r,2), (8) 

(32 = {h/2,)[{rs-l)/rs + tl/2] (9) 
MIa = 4/vIa = (1 + rl4)/rs, (10) 

where M2A is the outflow Mach number, V2a is the outflow Alfven speed, P2 = o,2s/'^2Ai 
and a2s is the outflow sound speed. Fig. 2 shows ti and M2A versus = P2/ Pi- Since 
t\ > 0, (8) shows that 2.5 < < 4 for a low Pi switch-off shock [7], with the lower limit 
being a perpendicular (±) shock and the upper limit a parallel (||) shock. As /3i — > 00, 

When B2 • V2/|-B2i'2| << 1, V2 will be supermagnetosonic [6] when = + > 
+ '^2A- Using (7), (2), (3), and (5) this condition reduces to Gr^ — 13rs — 20 < 0, and is 
satisfied for < 3.2 or ti > 1.25 from (8). A nearly uniform supermagnetosonic outflow 
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becomes the condition for a fast shock when the field is hne-tied at the outflow boundary. 
The jump conditions, (1) — (6), across such a quasi- _L fast shock for F = 5/3 give the 
equation 

M^A = 3r//?2/(4 - rf) + (3/2)r/(r/ - l)/(4 - r/). (11) 
Combining this with (8), (9), and (10) we obtain, 
rf = (12rs-30)-^[4r^-20rs + 7+(16r^-544r^+2952r^-4312rs-431+2400/rJ^/2]_ (^2) 

Fig 2. shows that 1 < r/ < 2 when 3.2 > > 2.5. The inverse dependence is expected 
because a decrease in corresponds to an increase in tension force along the shock plane, 
and thus a larger M^a, accounting for the larger r/. For the canonical quasi-± built-in 
slow shock with a line-tied outflow, a fast shock is likely, as the supermagnetosonic outflow 
condition requires only that ti > 1.25. 

The compression ranges of 2.5 < < 4 for the built-in slow shock and 1 < r/ < 2 
for the fast shock are important for shock acceleration theory: For distribution functions 
isotropic to first order in v/v* where v* is the particle velocity in the proper frame of the 
bulk flow V, the steady state Boltzmann equation can be written as a diffusion-convection 
(DC) equation. We deflne N{x,pp)dpp = 4:7rppf{x,pp)dpp, where / is the Boltzmann 
distribution function, x measures position, and Pp is the particle momentum. The DC 
equation across a general shock is then [9] 

dnKN - KndnN] - {l/3){dr,V^)dp^[ppN] = 0, (13) 

where is the normal flow velocity, is the normal diffusion coefficient, and we have 
assumed that gradients in the normal direction >> those along the shock. The solution 
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of (13) across the shock when the shock thickness << mean free path [9] shows that the 
outflow energy spectrum for a steeper inflow spectrum takes the power law form N oc 
Pp'^ with energy index w = (r + 2)/(r — 1) depending only on the compression ratio, r. 
Fermi acceleration operates as the particles diffuse between scattering centers (presumably 
turbulence) on each side of the shock. Particles always see the centers converging, as the 
normal velocity is larger upstream. 

Shock acceleration can dominate synchrotron loss when Tgyn, the shortest synchrotron 
loss timescale of the region, exceeds the longest shock acceleration timescale Tgh- 

Tsyn = 67rmeC/(7eSiCrT) > T^h ~ Kn\\/vl, (14) 

where ut is the Thomson cross section, 7e is the electron Lorentz factor, and is the 
diffusion coefficient normal to the slow shock and thus || to the downstream field. For 
particles moving at c, k^w ~ cA||/3, where A|| is the field gradient length [9] which we 
assume is of the same order in the inflow and outflow regions. From (14), the condition 
justifying the absence of a synchrotron loss term in (13) is then 

7e <~ .08{vi/ cm sec"^)2(Si/Gauss)"2(A||/cm)-^ (15) 

The third term in (13) can be thought of as the 1st order correction to the rest frame 
equation, when measured in the lab frame. This motivates the finding [10] that (13) 
includes not only guiding center diffusion through pitch angle scattering, but also drift 
acceleration from motion along the induced electric field. The per particle energy change 
from the latter increases with obliquity. For slow (fast) shocks, the curvature (gradient) 
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drift is 1 1 to the electric field and accounts for energy gains while the gradient (curvature) 
drift is anti-| I to the electric field and incurs particle energy losses [11]. These contributions 
conspire with those from the gyromotion component along the electric field for a net per 
particle momentum change dpp/dt = — PpVfj,, where v± is the fiow velocity ± to B. The 
relevant component of this force is included in the coefficient of the dN/dpp term of (13). 

The relative importance of a reconnection fast shock varies inversely with the size 
of the DR: When both shocks are present 3.2 > > 2.5 with 2.4 < < 3, so that 
1 < Tf < 2 with tu/ > 4, where Wg(^f^ is the slow (fast) shock energy index. (The range 
3.2 < Ts < 4 corresponds to a slow shock with no outfiow fast shock and 2.4 > Ws > 2.) 
Thus, if slow shock acceleration is effective, fast shocks cannot further steepen the already 
steep spectrum of particles that passed through the slow shock. However, if the DR length 
~ A||, then more of the fiow will see only the fast shock, and the spectrum from shock 
acceleration should have a somewhat lower energy index. 

Although we have used a DC equation, we recognize that shock acceleration is a rather 
non-linear process. However, fast shock simulations show that the Fermi acceleration 
engine is very efficient, transferring > 1/10 of the infiow energy to particles [9]. These 
particles tend to smooth out the shock by diffusion, produce turbulence, and increase the 
compression ratio above the jump condition value, as their escape and acceleration change 
the downstream equation of state. The shock smoothing can violate the assumptions 
built into the simple DC scheme, but the increase in the compression ratio over the linear 
limit enhances the non-thermal acceleration. The efficiency is relatively insensitive to the 
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obliquity of the inflow field unless the Mach number exceeds ~ 30. For our case, outflow 
fast shocks would then be effective if M2A ~ i^sti < 30 from (10). 

A recent hybrid simulation [12] of oblique slow shocks and the ion-ion cyclotron in- 
stability shows that steady or cyclically reforming slow shocks can develop depending on 
the inflow conditons. In the steady case a coherent Alfven wave train forms downstream. 
In the unsteady case the shock quasi-periodically transforms from a thin sharp transition 
to a wide diffuse transition, and Alfven turbulence is seen downstream. This cyclic refor- 
mation is strikingly similar to that seen in quasi- 1| fast shocks [9] and is an example of an 
electromagnetic beam instability brought on by the interaction of backstreaming ions with 
the inflowing plasma. Waves produced in the upstream by such instabilites are amplifled 
as they convect back to the shock front. The compressed waves then interact strongly with 
the inflow particles, scattering and slowing them, producing the entropy required for the 
shock. Some of the ions are scattered back upstream by the waves, and a subset of those 
are scattered back to the shock. The reformation of the thin structure occurs if the back- 
streaming particles leave the simulation region, and then no longer produce waves that are 
convected to the shock. The above process is in fact how non-thermal ions are extracted 
from an initially thermal input, initiating the Fermi process. It is the effectiveness of the 
Fermi acceleration which makes shock acceleration such a non-linear process. 

Slow shocks in the geomagnetic tail show wave substructures with properties similar 
to those of Ref. [12], and also show turbulence ahead and behind the shock fronts [13]. In 
addition, although much of the shock acceleration goes into ions, significant non-thermal 
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tails in the electron spectra are seen [14] , and are not modeled by hybrid simulations which 
assume a fluid electron population. Fast shock simulations show electron acceleration with 
injection electrons [15], and we expect that slow shocks could operate similarly. More 
simulations of slow shocks are needed which predict the spectrum of accelerated particles. 

Jet plasma in AG may be largely pair plasma [16], so ion-electron simulations might 
not be applicable. We suggest that reconnection and its shocks are a strong candidate for 
solving the re-acceleration problem in jets, which requires sustaining non-thermal electron 
emission over distances exceeding CTgyn, where c is the speed of light. Shock electron 
acceleration requires injection particles [16] which we now show, and argue that the DR 
may be able to provide such injection. 

We can estimate the range of particle energies accelerated by the Fermi process for the 
slow and fast shocks, and in particular the range of 7e. First consider the slow shock: An 
upper limit can be found by ensuring that the particles see an ordered held. This requires 
A|| > ^2 where g2 ~ c^7eme/(e-B2) is the particle gyroradius associated with the smaller 
field of the two flow regions. This implies 7e < eB2X\\/{meC^) ~ 10~^(S2/Gauss)(A||/cm). 
A lower limit can be found by demanding that the downstream particles be able to diffuse 
upstream. This requires particles of large enough energy to resonantly interact with the 
plasma waves which pitch-angle scatter the particles upstream. For Alfven turbulence 
[16], the electron lower bound is a factor ~ mp/mg times that for protons and is given by 
7e > {mp/me){v2A/c) where rup is the proton mass and e is the charge. Thus for the slow 
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shock 

{mp/me){v2A/c) < 7e < 10"^(S2/Gauss)(A||/cm). (16) 

Consider now the fast shock: Since the flow downstream from the fast shock is primarly 
_L to the shock normal, diffusion across the shock requires [9] 

Kn±/V3 > ^3 ~ 7emeC^/(eS3), (17) 

where vs is the downstream flow speed, is the downstream field, gs is the associated 
gyroradius, and is the diffusion coefficient normal to the shock but ~± B3. For A|| > ^3 
we use K,n\\ ~ cA||/3, so [9] = C5f|/(3A||) and (17) gives 7e > SeSsDsAn /(mgC^). This 
limit must be combined with the analogous limits as described for the slow shock so that 
for the fast shock we have 

Max [(mp/me)(f3A/c), 10~^(i?3/Gauss)(A||/cm)(^;3/c)] < 7e < 10~^(i?2/Gauss)(A||/cm), 

(18) 

where v^a is the Alfven speed downstream from the fast shock. Note that the upper limits 

in (16) and (18) are less than the upper limit in (15) when vf > (.01cm • sec~^Gauss~^)i?ii?2A|| .| 

The highly dissipative energy conversion in the DR could provide the injection elec- 
trons. To see this, note that upon absorbing the annihilated field energy, the average 7e 
there ~ (vf^/c^)(mp/2me). Combining this with (17) we see that slow shocks of large 
obliquity favor a DR which can inject, as the condition is 

via/v2A = r\l'^ jcx >~ cjvxA- (19) 
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We have discussed non-thermal acceleration by low f3i reconnection slow shocks and 
outflow fast shocks. The particle spectra across the slow shocks should be at least as flat as 
that given by the DC spectral index range of 2 < tUg < 3 for sufficiently large slow shocks. 
The dissipation region can provide injection electrons for acceleration and the above index 
range is consistent with observed features of radio galaxy lobes and jets [16]. 
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